Magnetic polarons in two-component hard core bosons 
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We use a high-temperature expansion to explore spin correlations around a single hole in a two- 
dimensional lattice filled with a hard-core two component bose gas. We find that the spins around 
the hole develop ferromagnetic order and quantify the degree of polarization at temperatures of 
order the hopping energy, finding a measurably nonzero polarization. We also discuss the effect of 
fixing the overall magnetization of the system for finite-sized systems. 
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In the mid 1960 's, Thouless and Nagaoka studied the 
two-component Fermi system on a bipartite lattice with 
very strong on-site repulsion [TJ [5] . They found that in 
the presence of a single hole the ground state was a fully 
polarized ferromagnet. These and further studies showed 
that at finite temperatures the system is not fully polar- 
ized: near the hole there is a ferromagnetic "bubble", 
while far away the spins are uncorrelated [3] |4j. On 
such bipartite lattices, the statistics are irrelevant for the 
single-hole problem, and the same physics should be seen 
in the bosonic case as in the fermionic system. Thus 
the bosonic ground state is the Nagaoka state, and at 
finite temperature, ferromagnetic correlations are found 
near the hole. Here we calculate these correlations in a 
two-component gas of hard core lattice bosons. We find 
that at experimentally relevant temperatures these cor- 
relations are measurable using a quantum gas microscope 

This is the simplest example of emergent physics in a 
strongly correlated system. Variants of it are also highly 
nontrivial: for example the ground state of two compo- 
nent fermions on a non-bipartite lattice with a single hole 
is unknown. A qualitative picture of this ferromagnetism 
can be developed by imagining a child's puzzle where tiles 
slide on a square grid. One tile is missing. By moving this 
"hole" one can rearrange the tiles. Here we have a quan- 
tum mechanical version of this puzzle. The motion of 
the hole from one location to another involves summing 
all possible paths. If the tiles are in a symmetric su- 
perposition of all possible arrangements (corresponding 
to ferromagnetism) then these paths will add construc- 
tively, allowing the hole to move over large distances. 
This ferromagnetic arrangement thereby minimizes the 
zero-point energy of the hole. 

Borrowing the term from how electronic motion cou- 
ples to lattice distortions, the elementary excitation con- 
sisting of a hole dressed by a ferromagnetic cloud is re- 
ferred to as a "polaron" . Other cold-atom polaron prob- 
lems include the behavior of a single down-spin atom in 
a Fermi sea of up-spins [SHU] . 

Even far from the strong-coupling hard-core limit stud- 
ied here the physics of two-component bosons is quite 
rich. This physics has been explored in theoretical works 



[j"UHT4"l , and in cold gas experiments [T51 US] ■ The com- 
ponents can be different hyperfine states [15], or differ- 
ent atomic species [15] • In the most ordered state there 
will be two independent order parameters, and it costs 
energy to twist the phase 4>i,(f>2, of either condensate. 
Depending on interaction parameters one can also find 
states where only some linear combination of the two 
phases has a finite stiffness. For example, with suffi- 
ciently strong attraction between the species there will 
be a condensate of "pairs" but no single particle conden- 
sate [T7]: One then has a stiffness to twisting cj>\ + <j) 2 , 
but not cf)\ — (f>2. Even more exotic is the "counter- 
superfluid" phase formed when the interspecies repul- 
sion becomes strong: One then has a stiffness to twist- 
ing </>i — 4>2- Under these circumstances trying to drive 
a current of species 1 to the right creates a current of 
species 2 to the left. Identifying the two components 
as the ±z-component of a pseudospin-1/2 object - the 
countcr-supcrfiuid state corresponds to an x — y ferro- 
magnet. If the in-species interactions are not sufficiently 
strong, either of these exotic states can be preempted by 
phase separation or collapse. [TS] In the single hole limit, 
the phase stiffnesses scale as the inverse of the system 
size. 

Here we use a high temperature expansion to calculate 
the correlations between spins bordering a single hole in 
a two-component hard-core Bose system on the square 
lattice. Using the techniques in [5j [T9H23] , these correla- 
tions can be directly measured, giving a signature of this 
interesting physics. 

The temperature scale at which these correlations be- 
come significant is of order the hopping energy t. In phys- 
ical units, this energy is on the order of t ~ ks X (1 nK) 
for 87 Rb atoms trapped by A = 820 nm lasers [53], but 
using lighter atoms such as 7 Li would increase the hop- 
ping energy energy and corresponding temperature by a 
factor of ten. Similarly, using a shorter wavelength lattice 
would also increase this scale. 

Our study assumes hard-core interactions, where dou- 
ble occupancy is forbidden. In most experiments, the 
strength of on-site interactions U is fixed and the hard- 
core regime is achieved by increasing the height of the po- 
tential barrier between neighboring sites so that t <C U . 
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Corrections to the hard-core results scale as t/U, Spiel- 
man et al. [24] report results with t/U ~ 0.001. 

Another relevant experimental detail is most cold atom 
systems are confined in harmonic traps. Local physics, 
such as the correlations we study, are unaffected by such 
confinement, as long as one restricts attention to regions 
where the polarons are dilute. 

The physics of Nagaoka ferromagnetism is relevant for 
a number of other cold atom systems [351 H3] . 

ANALYSIS 

We model the two-component Bose system via the sin- 
gle band Bosc-Hubbard Hamiltonian 

H = ~tJ2 E (4,<°crJ + a tj a <r,i) (1) 

(id) "=t4 

where a a ^ (aJ.J is the bosonic annihilation (creation) 
operator for a particles of type ("spin") a at lattice site 
i, and are all nearest-neighbor pairs and we limit 

ourselves to a two-dimensional square lattice. The single- 
particle spectrum has band-width of 8t. We work in the 
canonical ensemble, with fixed particle number, and do 
not need to include a chemical potential. 

The Bose-Hubbard model is a good description of the 
system as long as the band-spacing Ef, is large compared 
with the other relevant energy scales. We require t <C 
Eb so that the (single-particle) bands are distinct, while 
T -C Ef, is required so that all bosons are in the lowest 
band. In addition, we will be analyzing Eq. ([I]) within a 
high temperature expansion, requiring that the ratio T/t 
is not too small. 

For a cold-atoms experiment described by the single- 
band Hubbard model, the band spacing varies with mi- 
croscopic parameters as Eb ~ %/Vo-Er, where Vb is the 
height of the potential barriers between lattice sites, 
Er = h 2 k 2 /2m: k = 2tt/X being the laser wavenum- 
ber and m the particle mass. The tunneling t depends 
exponentially on Vb and is typically t ~ 0.1 — 0.01-Er for 
Er [27] , There is therefore a separation of scales, 
allowing t ~ T <C Eb- Deeper lattices accentuate this 
separation, at the cost of requiring lower temperatures. 

Strong interactions imply a hard-core constraint 

4,i4,< = °< ( 2 ) 

which is valid when the on-site interactions are large com- 
pared to t. 

We examine the case of a single hole in an infinite 
system and calculate the finite temperatures expectation 
values of an observable operator X by 

(X) = - Tr le^ 6 ; Z = Tr e~ fi6 (3) 



where = (ksT) is the inverse temperature; we take 
ks = 1. The trace is most readily calculated in a basis 
given by placing the hole on the site r^, and specifying the 
pseudospin Oi =j" / ^ on all remaining sites i ^ r^. We 
will look at the correlations between spins on positions 
which are fixed relative to ?v For observables of that 
form, denoting by £ a spin state with the hole at the 
origin, we have 

{X ) = N f^J2 X (0(Z\e-e 6 \0 (4) 
€ 

where the factor N S i te s comes from summation over all 
N S ites possible locations of r/i, and X |£) = X (£) 

To perform the calculation we use a high-temperature 

—j3H\ . Each power of 

H corresponds to a single "hop" of the hole, and the 
moments can be calculated from the sum of all closed 
paths of length n ( a n — paths" ) starting at the origin, 

^\( c PH) n \o = m n E (ein(f». ( 5 ) 

p£n— paths 

Here V v (£) is the spin permutation that results from 
moving the hole through the path p. Any open paths, 
that do not take the hole back to the origin, do not con- 
tribute to the sum, and the expectation value is zero if the 
path leads to a non-equivalent spin configuration. This 
requirement also restricts the sum to even values of n. 

Although the number of closed paths grows exponen- 
tially with n, we are able to exhaustively enumerate them 
for small n < 2M, and calculate a high temperature ap- 
proximant 

n=0 ^ 1 p£2n-paths £ 
n— pG2n— paths £ 

(6) 

We use M = 6. 

Estimating the error of cutting off such series to be on 
the order of the last term calculated, the correlation func- 
tions for spins around the hole are accurate to about 10% 
down to T/t ~ 0.4 for M = 6. To investigate lower tem- 
peratures, one would need to resort to more sophisticated 
methods of summing the series, such as the Monte-Carlo 
approach of Raghavan and Elser [55]. Lower tempera- 
tures are difficult to achieve experimentally. 

VACANCY-INDUCED FERROMAGNETISM 

The tendency towards ferromagnetism is apparent in 
the structure of Eq. . Ferromagnetic configurations £ 
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automatically have V (£) = £, regardless of the path p. 
A further insight is that it is only paths with loops in 
them that favor ferromagnetism. Paths p which retrace 
themselves have V (£) = £ regardless of £. 

To measure the polarization around the hole we define 



5's 



E S 



(7) 



where S\ is the spin operator applied to the site i and the 
summation is over the eight nearest-neighbor and next- 
nearest neighbor sites of the hole. The ground state of 
our system possesses a spontaneously broken symmetry 
In an infinite system with an infinitesimal magnetic field 
along z, Sg will have a finite expectation value. This ex- 
pectation value vanishes as T — > oo and approaches 4 as 
T — > 0. If there is no symmetry breaking field, then the 
spontaneous symmetry breaking occurs in a random di- 
rection. In a typical cold-atoms experiment, every time a 
new sample is created, this symmetry-breaking direction 
will be different. Under those circumstances, one can 
model the ensemble measurement by taking expectation 

values in zero field. By symmetry, in zero field (.Sg\ = 0, 
at all T, but the temperature dependence of its distribu- 
tion will be non-trivial. At T — > oo when all states are 
equally likely we expect a binomial distribution around 
zero. At T — ¥ 0, the distribution is uniform. This may 
be understood in several ways; in a quantum mechani- 
cal treatment, one would attribute this to the fact that 
each projection m of the spin multiplet is equally likely. 
Classically the z-component of a uniformly distributed 
random 3D unit vector is uniformly distributed. These 
distributions are shown in Fig. [T] 
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FIG. 1: The probability distribution of (5*8), the total spin of 
the bosons around the hole, at high and low temperature. 

To quantify these distributions, we examine the vari- 
ance of Sg. We define 



C« = 



Sg) -2 




E % 



i^n.n.n. 



i^n.n.n. 



(8) 



which is normalized and offset so that (Cg) goes to unity 
when the hole is maximally polarized and to zero when 
all sites are uncorrelated. Note that individual measure- 
ments of Cg can be negative or greater than one. 
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FIG. 2: (Cg,), the measure of polarization around the hole, 
and ACs = J (Cs 2 ) — (Cs), as a function of the relative tem- 
perature T/t. Note that ACs goes to ^9/28 as T — > oo. 



We have calculated for a range of temperatures (Cg 



and the uncertainty ACg = y (Cg 2 ) — (Cg) and they 
are shown in Fig. [2j In particular, at temperatures cor- 
responding to T/t — 0.4 we predict (Cg) = 0.05 and 
ACg = 0.62. This compares with a T — > oo result of 
(Cg) = and ACg = t/9/28 « 0.57. Both the non- 
zero mean of this quantity and the increase in variance 
are indicative of the ferromagnetic correlations present 
around the hole. About 5000 measurements would be 
needed to determine the mean to within 20% of the pre- 
dicted value. A given sample will contain multiple holes, 
so each experimental run can contribute multiple inde- 
pendent measurements. 



FIXED MAGNETIZATION 

In a cold atom experiment the number of f-spin and 
|-spin atoms are fixed, requiring a slightly different en- 
semble. This difference only matters when the correla- 
tion length becomes of the same order as the system size. 
For the temperatures described in Fig. [2] the correlation 
length is of order the lattice spacing, and these subtleties 
are irrelevant. 
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By using exact diagonalization on a small system we 
can, however, show that at an order of magnitude lower 
temperature one must consider these finite size effects. 
We consider a system of 5 x 3 sites described by the 
Hamiltonian in Eq. [TJwith periodic boundary conditions, 
7 f-spins, 7 |-spins and a single hole. We define a similar 
operator to the one used before 



f _ 
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56C,° 



(9) 



The constant C|° = (S^S^) is the infinite temperature 
two-spin correlation caused by the finite number of spins: 

^2° — I 2-n - t = f° r an ec L ua l number of f and J.- 
spins. 

The results are shown in Fig. [3} At high temperatures, 
one sees behavior indistinguishable from Fig. [2] while 
at low temperatures the expectation value is suppressed. 
This suppression can be attributed to the ferromagnetic 
order parameter being forced to lie in the x — y plane. 
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FIG. 3: (Color online) {C'i/ f° r a system of 3 x 5 sites with 
an equal number of f and 4- particles, as a function of temper- 
ature T/t. (solid blue line) Results from exact diagonaliza- 
tion, (dashed red line) result from high-temperature expan- 
sion taken to the same order as in Fig. [2] The two match well 
to about T/t ~ 0.4. 



OUTLOOK 

The problem of how charge and spin degrees of freedom 
interact with one another is key to a number of important 
condensed matter systems, most notably high tempera- 
ture superconductors. More importantly, conceptually 
clean examples of strongly correlated phenomena, such 
as the two component Bose system one, are essential to 
developing new paradigms for many-body physics. 

In a cold gas experiment the quantities (Ss) and (Cs} 
can be measured by a variant of the quantum gas mi- 
croscope technique pioneered by Bakr et al. [5j and ex- 
tended to spinor gases by Fukuhara et al 23J . An image 



is taken of the optical lattice, which shows the location of 
all particles, and their spin projection along a fixed axis. 
One would locate an isolated hole in this picture, and 
add up the spin projections of its neighbors to produce a 
single realization of Ss or Cs- The experiment would be 
repeated many times. A histogram similar to Fig. [T] can 
be produced for 5s- The ensemble average can be com- 
pared with our prediction for the quantum mechanical 
expectation value (Cs). 

While the single-hole problem studied here is already 
interesting, the many-hole problem is even more rich. At 
zero temperature, the system is both superfluid and ferro- 
magnetic. Kuklov et al. [I0HI2] have used Monte-Carlo 
methods to explore the relative strengths of superfluid 
and magnetic stiffnesses. Although no finite tempera- 
ture studies have been done, both orders will disappear 
as one heats the system. It would be interesting to know 
if magnetism or superfluidity vanish first, or if the two 
orders vanish simultaneously [29]. This question could 
be largely answered by studying the interaction between 
two polarons. 
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